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To establish the bond-site duality of explosive percolations in 2 dimension, the site and bond 
explosive percolation models are carefully defined on a square lattice. By studying the cluster 
distribution function and the behavior of the second largest cluster, it is shown that the duality in 
which the transition is discontinuous exists for the pairs of the site model and the corresponding bond 
model which relatively enhances the intra-bond occupation. In contrast the intra-bond-suppressed 
models which have no corresponding site models undergo the continuous transition and satisfy the 
normal scaling ansatz as ordinary percolation. 
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Recently Achlioptas process (AP) [l[ which was 
suggested to show a supposedly first order transition on 
the complete graph triggered intensive studies on explo- 
sive percolations |2|-l7| . However subsequent studies have 
proved that the transition in AP on the complete graph 
is continuous 0-0- We have also shown that AP on the 
Bethe lattice shows a continuous transition Q . Therefore 
the transition in the original AP is physically established 
to be continuous in the mean-field level or in high dimen- 
sions. 

Until now studies on AP have been done mainly on the 
complete graph. Even though there are some studies on 
AP in 2 dimension (2d) 0, SS-Q, 

the transition na- 
ture in lower-dimensions is still not fully understood. For 
example, the bond percolation under AP with a product 
rule was first argued to show a discontinuous transition 
[H, However, based on the measurement of the largest 
cluster distribution, Grassberger et al. 0] argued that 
the bond percolation with the same product rule in 2d 
undergoes continuous transition 0]- The site percola- 
tion under AP with a product rule in 2c? has been first 
proved to undergo the discontinuous transition based on 
the detailed analysis of cluster size distribution and hys- 
In contrast Bastas et al. argued that the 
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site percolation under AP with a sum rule in 2d under- 
goes continuous transition based on the finite-size scaling 
analysis with relatively small system sizes fl2| . 

Such controversies 0, 0, 0, EH, E3 also imply that there 
doesn't seem to exist the bond-site duality among explo- 
sive percolation models on 2d lattices unlike ordinary per- 
colation jl3 , 14 1 . Here the bond-site duality means that 
a bond percolation model has the same transition nature 
or belongs to the same universality class as the corre- 
sponding site percolation model except for the properties 
depending on details of models such as transition proba- 
bility p c , etc. Moreover, if the transition is truly discon- 
tinuous, then determination of critical exponents from 
finite size scaling analysis like in Ref. [12j has no physi- 
cal meaning as we already addressed in Ref. [ll|. Thus, 
resolving such controversies by establishing the bond-site 



duality in 2d is theoretically very important and interest- 
ing. 

The controversies should come from ambiguities in the 
definitions of explosive percolation models on lattices. 
Therefore it is very important to make clear definitions 
of models with various growth rules for AP. There can 
be six kinds of models on a 2c? square lattice. Among 
them, two pairs are bond models. One pair consists of 
the bond models which physically enhance occupation of 
intra-cluster bonds. The other pair consists of bond mod- 
els which relatively suppress occupation of intra-bonds. 

In site percolation there cannot be the distinction be- 
tween inter-sites and intra-sites to a cluster. So there is 
no ambiguity in the definition of site models as in bond 
models. As we shall see, the intra-bond-enhanced mod- 
els and the corresponding site models show the bond-site 
duality in which the transition is discontinuous. In con- 
trast the intra-bond-suppressed models show the contin- 
uous transition. Physically there should be no site model 
corresponding to such intra-bond-suppressed models. 

There are two fundamental percolation models on 
lattices 13[ . One is the site percolation model and the 
other is the bond percolation model. In the site perco- 
lation, there is no new site occupation which does not 
change the size of clusters. Under AP, two vacant sites 
A and B are randomly selected. Let {s^} (or {sbj}) 
be the sizes of tia (or hb) clusters which would be con- 
nected by occupying the site A (or B). In the square 
lattice tia (or ub) is at most 4. In the Site model 
with a Product rule (SP model) the site A is oc- 
cupied if 1 x rj™=i s a, < 1 x n"=i s s. ■ Otherwise, 
the site B is occupied. Similarly, in the Site model 
with a Sum rule (SS model) the site A is occupied if 

(i+Er= A i^j<(i+E-ii^)- 

We define four bond percolation models under AP. 
First two unoccupied bonds a and b are selected ran- 
domly. If the bond a is an inter-bond, then it connects 
two different clusters of sizes s a i and s a 2- Then under 
a product rule the product £ a for the bond a is clearly 
defined as £ a = s a i x s a 2 without any ambiguity. If the 
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bond is an intra-bond, it internally connects two sites in 
the same cluster. Then £ a can be defined in two differ- 
ent ways. One is £ a = s al x 1 (Bond Product Type 
1 model: BP1 model). The other is £ a = s al x s al 
(Bond Product Type 2 model: BP2 model). The 
product £{, for bond b is similarly defined. Then occupy 
bond a if £ a < £&. Otherwise occupy bond b. Therefore 
two bond product models, BP1 and BP2, come from the 
ambiguity to define the product for the selected intra- 
bond. Similarly we can define two kinds of bond models 
with a sum rule. The sum a a for the inter-bond a is de- 
fined clearly as a a = s a i + s a 2- In contrast <j a for the 
intra-bond a can also be defined in two-different ways: 
a a = s a i+0 (Bond Sum Type 1 model: BS1 model) 
or a a = s i + s a i (Bond Sum Type 2 model: BS2 
model). Then occupy bond a if a a < a . 

The physical meaning of models is that type 1 models 
(BP1 and BS1) relatively enhance the intra-bond occupa- 
tion, whereas type 2 models (BP2 and BS2) suppress the 
intra-bond occupation. Thus if there exists the bond-site 
duality, it should be between type 1 bond models and 
site models. As we shall see, the duality exists for the 
pair of BP1 and SP models and the pair of BS1 and SS 
models. BP2 and BS2 have no corresponding site models 
for the duality. 

Until now, only the BP2 model has been studied for 
explosive bond percolation model 0, 0, 0] ■ In Ref. [§] 
BP2 model was argued to show discontinuous transition, 
whereas in Ref. [7( the same model was argued to un- 
dergo continuous transition. BP1, BS1 and BS2 models 
have never been studied until now. Both SP and SS 
have been studied as explosive site percolation models. 
SP model ll| has been proved to show discontinuous 
transition, whereas SS model [13] was argued to show 
continuous transition based on the numerical studies of 
relatively small systems. 

To understand the transition nature of percolation 
physically, the cluster size distribution should be the first 
one to understand. The cluster size distribution P s (p) at 
an occupation probability p of a bond (or site) is the 
probability that an occupied site belongs to a cluster 
which has s sites. It has been shown that P s {p) pro- 
vides an excellent method to determine p c as well as the 
transition nature for ordinary percolation (OP)[l3} 
and explosive percolations 0, [ll|. Thus, we first investi- 
gate P s (p) for each model to obtain p c and the transition 
nature physically. When p < p c , P s for OP decays expo- 
nentially as s increases. This means that the probability 
to find a large cluster vanishes exponentially. On the 
other hand, when p > p Cl there exists the macroscopi- 
cally large cluster (LC) and P s for finite s also decays 
exponentially with a peak for LC. At p — p c it is well 
known that P s satisfies a power-law, P s (p c ) ~ s~ s , with 
5 ~ 1.055 for OP. 

In contrast, P s for SP model in Fig. HJa) is completely 
different from that for OP as shown in Ref. [ill. When 
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FIG. 1. (Color online) Plots of P a (p) against s for various 
models. Unfilled symbols denotes P a for the largest cluster 
on the lattice with the linear size L. (a) SP model at p — 
0.7723(7) (w p c ). Insets: at p = 0.75 (< p c ) (upper inset) 
and p = 0.78 (> p c ). The line with S = 0.95 means the 
relation P s ~ s -0 ' 95 . (b) BP1 model at p = 0.6937(8) (« p c ). 
Insets: at p = 0.67 (< p c ) (upper inset) and p = 0.70 (> p c )- 
The line denotes the same thing as in (a), (c) BS1 model at 
p = 0.5979(4) (« p c ). Insets: at p = 0.590 (< p c ) (upper 
inset) and p = 0.605 (> p c ). The line denotes the relation 
P s ~ s^ 1 - 5 (d) SS model at p = 0.6916(5) (» p ). Insets: 
P 3 at p — 0.685 (< p c ) (upper inset) and p = 0.695 (> p c )- 
The line denotes the same thing as in (c). (e) BP2 model at 
p = 0.5266(1) (» p c ). Insets: at p = 0.5250 (< p c ) (upper 
inset) and p = 0.5270 (> p c ). The line means the relation 
P s ~ s" 102 . (f) BS2 model at p = 0.5270(1) (« p c ). Insets: 
at p = 0.5250 (< p c ) (upper inset) and p = 0.5275 (> p c )- 
The line denotes the same thing as in (e). (g) Scaling collapse 
for the relation {TJ of BP2 model and (h) the collapse of BS2 
model. Used exponents in (g) and (h) are 5 = 1.02 and vo = 
0.51. 



p < p c , P s for SP model has a stable hump in the tail as 
p — i- p~ 0, [ll[ . The behavior of hump for SP model have 
been studied in detail in Ref. [ll[. On the other hand, at 
p ~ p c P s for LC starts to be detached from the continu- 
ous distribution of P s for finite s, which satisfies P s ~ s~ s 
with S = 0.95(1)(< 1). Such power-law behavior for fi- 
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nitc s is observed for sufficiently large p(> p c ) (see the 
lower inset of Fig. [T{a)) Based on this typical be- 

havior of P s (p) we determined p c for SP model [ljj . Since 
the hump contribution for p < p c does not depend on the 
lattice size L and S < 1 at p c , there should be many sta- 
ble large (but still microscopic) clusters before transition, 
which stro ngly indicates the discontinuous transition for 
SP model [ll]. To determine the p c 's and find the tran- 
sition nature for other models, we now analyze P s 's as in 
Ref. [HI. 

In Fig. mb) we display P s (p) for BP1 model. The 
data for p < p c in the upper inset of Fig. QJb) clearly 
shows that the hump in P s (p < p c ) docs not depend on 
L as for SP model in Fig. [TJa). At p ~ p c , P s for LC 
starts to be detached from the continuous distribution of 
P s as shown in the main plot of Fig. QJb). From the 
best fit of the data for finite s, we obtain P s <~ s with 
5 ~ 0.95(1)(< 1). For p > p c we find that P s for finite 
s still satisfies the power-law with 5 < 1 (see the lower 
inset). These behaviors of P s (p) for BP1 model exactly 
coincides with that for SP model, which is known to un- 
dergo discontinuous transition [Uj. The only difference 
between BP1 model and SP model is in the value of p c : 
p c ~ 0.6937(8) for BP1 model and p c ~ 0.7723(7) for 
SP model. The coincidence of P s between BP1 model 
and SP model physically means that there is the bond- 
site duality under AP with the product rule if the bond 
model enhances the intra-bond occupation. 

Figs. \Hc) and (d) show P s for BS1 and SS models. By 
assuming the relation P s ~ s s only for small s at p ~ p c , 
we obtain 5 ~ 1.5 for BS1 and SS models. However, as 
depicted in Figs. QJc) an d (d), P s for BS1 and SS models 
seems to substantially deviates from the power-law due to 
a main contribution from the hump distribution for large 
s. Even for p c at which P s for LC starts to be detached 
from the seemingly power-law like regime, the hump still 
exists. The contribution of the hump distribution part to 

L 2 L 2 

J2 S P{ s ) = P f° r site percolation (or J2 s =i P( s ) = 1 f° r 
bond percolation) even at p = p c is more than 90%. This 
means that the contribution of the power-law like regime 
for small clusters to ^ s P(s) is trivial. This result im- 
plies the absence of the singular behavior, P s (p) ~ s~ s 
for s — > oo at p c , which is the intrinsic property of the 
continuous transition. Moreover in both BS and SS mod- 
els the hump distribution for p < p c does not depends on 
the lattice size L (see the upper insets). This anomalous 
behavior suggests that there exists a type of powder-keg 
1 51 ] in BS1 and SS models, unlike in the fully connected 
networks Q. In addition, the existence of the hump dis- 
tribution even for p > p c (see the lower insets) implies 
the existence of multiple stable macroscopic clusters for 
p > p c in BS1 and SS models. Furthermore, P s for BS1 
model (Fig. QJc)) is nearly identical to P s for SS model 
(Fig. Old)) as the case for BP1 and SP models. The 
coincidence of P s between BS1 and SS models except for 



the value of p c also implies that there is the bond-site 
duality under AP with the sum rule if the bond model 
enhances the intra-bond occupation. 

In Figs. [He) and (f) we display P s for BP2 and BS2 
models. P s for BP2 and BS2 models is physically differ- 
ent from P s for BP1 and BS1 models. Even though P s 
for BP2 and BS2 models has a hump in the tail when 
p < Pc (upper inset), this hump vanishes in the limit 
L — > oo and P s at p c shows the power-law singularity 
Ps{p) ~ s~ s with S = 1.02(1) > 1. As we increase p 
further to p > p c , P s for finite s decays exponentially as 
in OP (see the lower insets for p > p c ). This P s behav- 
ior is physically the same as that in OP. As for OP, P s 
or n s (= P s / s) at p c for the percolation with continuous 
transition satisfies the scaling relation [l3[ 

P s = s- 5 f(s/L l /" v ) (n s = s-^fis/L 1 ^)), (1) 

where 5 = t — 1 and v is the correlation length critical 
exponent. As shown in Figs, [2(g) and (h), P s (or n s ) for 
both BP2 and BS2 models satisfies the scaling relation 
(HI) with S = 1.02(t = 2.02) and va = 0.51 very well. P s 
in Figs. H](e), (f), (g) and (h) strongly suggests that both 
BS2 and BP2 models show continuous transition. Our 
estimation of p c = 0.5266(1) for BP2 model is nearly 
the same as those in Refs.@, Q, even though the two 
references were contradictory to each other in transition 
nature. 

From the obtained r and a from Eq. ([T]), one can 
calculate the critical exponents (3, 7 and is, which must 
be identical to the values evaluated from the finite-size 
scaling (FSS) properties of the average size S(p,L) of 
the finite clusters and the order parameter Poo(p,L) if 
the transition is truly continuous. From S(p, L) for BS2 
model obtained by the numerical simulations with L = 
512 — 4096 and the relation p m ax(L) = p c + bL~ 1 / 1 ' , where 
Pmax(L) is p at which S(p, L) is maximal [l3j |. we obtain 
v = 1.00(1) and p c = 0.5270(1) for BS2 model. From the 
relation for the maximal value of S(p,L) : S max ~ L 1 ^ 
[l3| and the obtained data for S we obtain 7/1/ = 1.90(2) 
Thus S{p,L) for BS2 model satisfies FSS ansatz [13j 
S(p, L) = Li/vfdp-pjL 1 ^) very well with v = 1.00(1) 
and p c = 0. 5270(1). Poo (p, L) for BS2 model satisfies the 
FSS ansatz ^ Pao(p,L) = L-^ u g((p - p^L 1 /") very 
well with (3/v = 0.04(1). Since 7 and (3 are related to r 
and a as 7 = (3 — t)/o and /3 = (t — 2)/<j [l3| , v, 7 and /3 
13( obtained from FSS ansatz are consistent with r and 
a obtained from Eq. (fTJ for BS2 model. The obtained 
exponents are r = 2.02(1), a = 0.51(1), v = 1.00(1), 7 = 
1.90(2) and /? = 0.04(1). For BP2 model we obtain iden- 
tical exponents and scaling relations to those for BS2 
model. 

The discontinuous transition is characterized by the 
existence of metastable states. The metastable states in 
percolation transition are generally originated from the 
coexistence of multiple macroscopic clusters. To study 
multiple macroscopic clusters, we now focus P2(p, L) with 
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FIG. 2. (Color online) (a) Plot of P2{p, L) against p (a) for or- 
dinary bond percolation, (b) for SP model(Inset:BPl model), 
(c) for SS model (Inset:BSl model) and (d) for BS2 model 
(Inset: BP2 model). 



which a randomly chosen occupied site belongs to the 
second largest cluster. If the transition is discontinu- 
ous, then P2 approaches a nonzero value at p ~ p c in 
the limit L — >• 00. On the other hand, when the transi- 
tion is continuous, the largest cluster grows by gradual 
adding of small clusters and P2 at p ~ p c should de- 
crease to zero in the thermodynamic limit. For the sake 
of comparison, we first measure Pi for ordinary bond 
percolation (OBP). As shown in Fig. GJa), p at which 
P2 is maximal, P2max, approaches to the known value of 
p c = 1/2 as L increases. Furthermore the maximal value 
Piipimax) decreases as L increases. These behaviors of 
P2{p, L) clearly show the absence of the multiple macro- 
scopic clusters in OBP for p ~ p c in the limit L — > 00 
and the transition becomes continuous. 

In Fig. fflb) we show P 2 (p, L) for SP and BP1 models. 
Here we also see the bond-site duality for BP1 and SP 
models. Like P2 in OBP, P2max for BP1 and SP models 
approaches to the estimated p c from Figs. Ufa) and (b) 
as L increases. However, in contrast to OBP, P2(p2max) 
for both models does not decrease and remains nearly 
at constant value as L increases. This behavior clearly 
shows that P2 for both models does not vanish at p ~ p c 
in the limit L — » 00, which provides a strong evidence for 
the discontinuous transition. 

In Fig. [2jc) we display P%(jp, L) for BS1 and SS models. 
As can be seen from Fig. Efc), there also exists the bond- 
site duality for BS1 and SS models. As shown in Fig. 
[He), P2 for both BS1 and SS models manifests anomalous 
behavior. Unlike OBP, p^max for both models hardly 
varies as L increases and is very close to the estimated 
p c from Figs. QJc) and (d) regardless of L. Furthermore 
P%{p2max) for SS model remains nearly at constant value 
or increases slightly as L increases. These results for 
BS1 and SS models physically mean that there exists a 



stable macroscopic second largest cluster even at p c in 
the thermodynamic limit and the transition should be 
discontinuous. 

In contrast, P 2 for BP2 and BS2 models in Fig. gfd) 
is physically very similar to that for OBP in Fig. Ufa). 
Vimax for BP2 and BS2 models approaches to the esti- 
mated pc as L increases. P 2 {p2max) for BS2 model de- 
creases as L increases as that for OBP, which indicates 
that P 2 vanishes at p ~ p c in the thermodynamic limit. 
Thus, the transition becomes continuous for both BP2 
and BS2 as expected from P s (p) data. 

In this letter we exactly define the explosive lattice 
percolation models on the square lattice. By studying P s 
and P 2 for the models, we observed the bond-site duality 
in the pair of SP and BP1 models and in the pair of SS 
and BS1 models. The duality means the discontinuous 
transition. In contrast two bond models, BP2 and BS2 
models, which relatively suppress the intra-bond occupa- 
tion, undergo the continuous transition, which satisfies 
the normal scaling behavior like Eq. (1). 
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